Interplay of quantum and thermal fluctuations in a frustrated magnet 



o 
o 

(N 
D 

in 



S. V. Isakov^ and R. Moessner^ 
^Department of Physics, Stockholm University, AlbaNova, S-106 91 Stockholm, Sweden and 
^Laboratoire de Physique Theorique de I'Ecole Normale Superieure, CNRS-UMR8549, Paris, France 

We demonstrate the presence of an extended critical phase in the transverse field Ising magnet 
on the triangular lattice, in a regime where both thermal and quantum fluctuations are important. 
We map out a complete phase diagram by means of quantum Monte Carlo simulations, and find 
that the critical phase is the result of thermal fluctuations destabilising an order established by the 
quantum fluctuations, ft is separated by two Kosterlitz-Thouless transitions from the paramagnet 
on one hand and the quantum-fluctuation driven three-sublattice ordered phase on the other. Our 
work provides further evidence that the zero temperature quantum phase transition is in the 3d XY 
universality class. 
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I. INTRODUCTION 

In geometrically frustrated magnets, the arrangement 
of the spins on a lattice precludes the establishment of 
a simple coUinear Neel state. In the absence of such an 
'obvious' ordering pattern, these magnets are often sub- 
ject to a large number of competing instabilities. This 
endows them with a rich behaviour, manifested in un- 
conventional phases, phase transitions and excitations, a 
large low-energy density of states and unusual dynamical 
properties^ 

When added to a frustrated magnet which is not or- 
dered classically at temperature T = 0, both thermal and 
quantum fluctuations often generate an ordering transi- 
tion, a process known as order by disorder. Whereas 
quantum order by disorder occurs in a fashion quite anal- 
ogous to the ordering induced by thermal fluctuations, 
the nature of the ordering (if any) of even some 'simple' 
model systems (such as the Heisenberg magnet on the 
pyrochlore lattice) in the presence of stronger quantum 
fluctuations remains unclear. 

The interplay of thermal and quantum fluctuations has 
so far not received a great deal of theoretical attention 
for frustrated magnets. By contrast, there has been a 
considerable amount of interest devoted to the real-time 
dynamics and transport at finite temperature near quan- 
tum critical points.^ In addition, this combination has 
even been claimed to provide an alternative use of quan- 
tum effects for speeding up computations!^ 

In this publication, we study the combined effect of 
thermal and quantum fluctuations on possibly the sim- 
plest quantum frustrated magnet, namely an Ising model 
on the triangular lattice in a transverse magnetic field. 
The manifold of classical ground states is hugely degen- 
erate, and correlations averaged over it are criticalJ This 
model thus enables us to investigate how the two types 
of fluctuations together establish and destroy order out 
of the exponentially large set of classical ground states. 
By changing the values of transverse field, F, and tem- 
perature, T, their strengths are in principle separately 
experimentally tunable. 

We map out the phase diagram of the Ising model on 
the triangular lattice in the transverse field-temperature 



plane. The bulk of the paper is devoted to a numeri- 
cal study based on a continuous time quantum Monte 
Carlo algorithm. We find that there are three different 
phases. In the regime where quantum fluctuations are 
weak, yet dominate over thermal fluctuations, they gen- 
erate a three-sublattice ordered phase. Upon increasing 
the strength of thermal fluctuations, it 'melts' into a crit- 
ical phase which has a finite extent and displays drifting 
exponents. This dome-shaped critical phase is delineated 
above and below by Kosterlitz-Thouless (KT) transitions 
at finite temperature. At T = 0, it terminates in the clas- 
sical critical point at zero field and in a quantum phase 
transition at a finite Tc- At high temperatures and large 
field strengths, one retrieves a simple paramagnet - or- 
der by disorder is a reentrant phenomenon. We briefly 
comment on the advantages of different diagnostics for 
the floating KT phase. 

The results reported here largely build on the follow- 
ing pieces of work. In Ref. |^ a number of frustrated 
transverse field Ising models were studied, and a connec- 
tion was pointed out between these models and classical 
stacked magnets in a scaling limit. An influential pa- 
per by Blankschtein and coworkers^ on stacked triangu- 
lar magnets derived a Landau- Ginzburg- Wilson theory 
for these systems, and predicted the stacked problem to 
be in the universality class of the XY model in three di- 
mensions with a six-state clock term, a problem which 
has in turn been studied by Jose et al}'^ Collecting these 
results together led to a conjecture of the phase diagram 
we map out in the following We note that, by adding 
a continuous degree of freedom to a classical triangular 
Ising antiferromagnet, a floating critical phase has previ- 
ously been obtained by Chandra et ali^ 

The general validity of the Landau- Ginzburg approach 
has been somewhat called into question by Monte Carlo 
simulations on the ferromagnetically stacked triangu- 
lar Ising antiferromagnet, with several groups present- 
ing evidence in favour of critical behaviour in the three- 
dimensional universality class, while others found a dif- 
ferent, possibly new, universality class to be a more likely 
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Our work here lends further support 
to the Landau-Ginzburg approach which is independent 
of the previous numerical simulations. 
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II. MODEL 

Possibly the simplest model which has both a high clas- 
sical ground state degeneracy and a non-trivial quantum 
dynamics is the transverse magnetic field antiferromag- 
netic Ising model on the triangular lattice. The model 
has the following Hamiltonian: 

= J^affT|-r^af, (1) 

where J > is nearest-neighbour exchange couplings and 
r is the magnetic field strength. 

The classical model triangular Ising magnet is obtained 
in the absence of a transverse field, F = 0. Its solution is 
well knownii it is disordered at any finite temperature, 
and critical at T = 0, where it retains an extensive zero 
point entropy, S, with S/N = 0.323^3, where N is the 
number of spins and ks denotes the Boltzmann constant. 
Upon application of an infinitesimal transverse field, F — 
0+, at T = 0, the magnet orders into a three sublattice 
pattern with a sublattice magnetisation pattern 
(1,0,-1)^ 



A. Mapping to a stacked magnet 

Using the Suzuki- Trotter formalism^LiS one can map 
the d-dimensional quantum model onto the (d+1)- 
dimcnsional classical problem. The partition function 
of the quantum Hamiltonian (Eq. Trexp(— is 
equivalent to the following partition function: 

Z = Trexp(-S'c;), 

Scl = ^'^{iJ).T Si.rSj^T ^ X/i,T Si.TSi,r+l, (2) 

where Si^r are classical Ising spins taking values 1 or — 1; 
T is the index in the imaginary time direction and runs 
from to iV, N = /3/e is the number of slices in the 
imaginary time direction and e is the imaginary time 
step; K = Je and K'^ = — ^IntanheF. The mapping 
becomes exact in the scaling limit e — > 0, K'^ —^ oo, while 
maintaining 

e-^^^ = eF . (3) 

It is in this formulation that the dimensional crossover 
between the high temperature d = 2 and the low tem- 
perature c? = 3 behaviour is most transparent. Besides 
providing an efficient Monte Carlo algorithm (see be- 
low), the above mapping also allows us to make contact 
with the literature on stacked triangular magnetsjiS from 
which our problem can be obtained in the scaling anal- 
ysis of infinite coupling strength in the third dimension, 
with the third dimension being of finite length at nonzero 
temperature. 



B. LGW analysis 

Blankschtein et almi have identified a complex XY 
order parameter for a stacked triangular antiferromag- 
net: the amplitude and phase of the Fourier mode at 
Q± = (±47r/3,0), Vi = mexp(±i(^). They have ob- 
tained the following LGW Hamiltonian (see below): 

+U(im^ + vem^ cos(66'). (4) 

A six-fold clock anisotropy thus appears at sixth order 
in m. The six- fold clock term is irrelevant at d = 3 so 
that the transition into a phase with nonzero m should 
be in the d = 3 XY universality class. This applies both 
to the zero temperature quantum phase transition as we 
vary F and to the corresponding transition in a stacked 
Ising magnet. With the clock term being dangerously 
irrelevant, the transition is immediately into a state with 
a broken clock symmetry, the details of which depend on 
the sign of vq (see below). 

At finite temperatures the triangular transverse field 
Ising model maps onto a stacked magnet of finite size, 
I/T, in the temporal direction, which is effectively 2- 
dimensional as the temporal correlation length, ^r, ex- 
ceeds Lt- close to a continuous phase transition. The 
2-dimensional six-state clock model, studied by Jose et 
aZ.,^° has a remarkable phase diagram consisting of three 
phases: a disordered phase, an extended intermediate KT 
phase, and an ordered phase. The extended KT phase 
owes is existence to the fact that in d = 2, the clock 
term does not become relevant immediately below the 
KT transition of the XY model; instead, the coupling 
needs to be increased further until the clock term is suf- 
ficiently strong to assert itself. 

One of the central objectives of this paper is to estab- 
lish the presence of this floating KT phase as a finite- 
temperature induced dimensional crossover in the 2-1-1 
dimensional quantum Ising magnet on the triangular lat- 
tice. In our context, the presence of such a phase would 
be the result of an XY order parameter generated by frus- 
tration and stabilised by quantum fluctuations, together 
with destabilising fluctuations provided by the finite tem- 
perature. 

Details of the ordered state depend on the sign of 
vq. Hlcw is minimised by M = |M|exp(j$), with 
$ = rnr/S and <i> = {n + l/2)7r/3 for vq < and 
f 6 > respectively, with n = 1 ... 6. The real space 
configurations of the ordered phases can be obtained by 
Fourier transforming these modes. In the case of vq < 0, 
Fourier transforming yields the following magnetisations 
assigned to the three sublattices of the triangular lattice: 

|A/|(1, -1/2, -1/2). We depict this phase as (H ). In 

the case of wg > 0, we have the following sublattice mag- 
netisations: |M|(y3/2,0, -\/3/2). We depict this phase 
as (-1-0—). There are six degenerate states for both < 
and ug > 0. 
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It is difficult to determine the sign of vq - in par- 
ticular, in the effective Hamiltonian, its value may 
drift and even change signf2*i^ The partition function 
Z — exp{—J2ijSif^ijSj) of the Ising model H = 
Si^ij-^j is equivalent to the following partition 
function (see, e.g. Ref. l20[l 

Z oc y Dtp exp {— C) = J Dtp expl—'^^^iKijipj 

ij 

+ ^lncosh(^2K,,V,)}, (5) 



The first term is the interaction matrix, which de- 
termines the soft modes ip± to be located at Q± 
(Eq. [TJ. Performing the Fourier transform and expand- 
ing Incosh(x) — ^ — J2 l5 ^ Mfo + • ■ -I we obtain the 
following expression for the contribution of the non- linear 
term to C 



Cnl = -2^i^(k)if(-k)V(k)V^(-k) 
" 4 

[](if(k0^(k,))5'(ki + ... + k4) 

=1 

" 6 

^(^(k^)V'(k^))(5'(kl + ...-|-k6 



IE 



64 ^ 
45 ^ 



(6) 



ki U=l 



The interactions for the stacked triangular antiferro- 
magnet are written in Fourier representation as 

K(k) = (^J cos{k^)+cos{k^/2 + V3ky/2) 

+ cos{K/2 ~ Viky/2) - J' cos{k^)^ , (7) 

where the sum is over the hexagonal Brillouin zone. The 
LGW Hamiltonian is thus constructed in terms of the 
coefficients m exp(±i(/))), varying slowly in space, of the 
'soft' modes i/;±. 

The primed delta functions indicate that the wavevec- 
tors ki need to add up to a reciprocal lattice vector. This 
requirement leads to the XY nature of the effective theory 
at low order. At sixth order, ^'+'0- a-nd + occur, 
the latter being the clock term, mcos(66'). < at this 
order. However, the cos(60) term appears at higher or- 
ders as well, e.g., ■0+ '(/'-(■!/'+ + 0-) eighth order. Since 
the terms in the series have alternating sign and large 
coefficients, one cannot reliably determine the sign of ug 
near the transition this way. 



problem'. It can therefore be studied reasonably straight- 
forwardly by Monte Carlo simulations. We thus simulate 
the classical problem defined by Eq.|21 However, it is dif- 
ficult to simulate the discretized version of Eq.|21because 
of the scaling limit (Eq. O : in order to avoid discreti- 
sation errors, one has to take a very large ferromagnetic 
coupling in the imaginary time direction K"^ ^ oo at the 
same time increasing the height of the system exponen- 
tially as e^^^ . 

In order to avoid this problem, we use a continuous 
time algorithm. ^-'^•^^■^'^ The basic idea behind this algo- 
rithm is that, in the scaling limit, the density of domain 
walls in the imaginary time direction becomes exponen- 
tially sparse and it is thus more efficient to keep track 
of the location of the domain walls, using exp(2ii"^) 
as a unit of length. Thus, the height of the system, 
iV/exp(2iv:'') = l3/{eexp{2K^)) = /3r, measured in units 
of exp(2i4r '^), remains fixed in the continuum limit. 

Due to the frustrated nature of our problem we can- 
not use a cluster algorithmSli2Si2& in the space directions 
as one could for the case of unfrustrated magnets; how- 
ever, we can use a cluster algorithm in the time direc- 
tion. The algorithm works as follows (for more details, 
see Ref. |23) . We pick a random site on the triangular 
Lx L lattice and build a cluster on that site in the imag- 
inary time direction. Its length, t, is distributed accord- 
ing to the probability distribution P(t) oc exp(— Ft) - 
this prescription eliminates a 'freezing' of the algorithm 
due to the diverging temporal coupling. A given clus- 
ter is fiipped using the Metropolis prescription in the 
spatial direction, i.e. we fiip the cluster with probabil- 
ity p = min(l, exp(— Ai?)), where AE is the (spatial) 
energy difference between the original configuration and 
the configuration with a flipped cluster. 

One check of our Monte Carlo algorithm consisted of 
comparing it with the diagonalisation of a 3 x 3 lattice, 
and we have found excellent agreement between the two 
approaches. 



B. Parameters of the simulations 

We impose periodic boundary condition on the trian- 
gular lattice and performed simulations on lattices of size 
i = 9, 12, 15, 18, 24, 30, 36, 48, 60, 75. We have estimated 
the correlation time Tc for different parameters. We per- 
form usually 500tc Monte Carlo steps (MCS) for equili- 
bration and from IO^'tc to 5 • IO^Tc MCS for averaging. 



III. MONTE CARLO METHOD 

A. The quantum Monte Carlo algorithm 

The cZ = 2 -I- 1 dimensional classical problem obtained 
this way has only positive weights - there is no 'sign 



IV. SCALING ANALYSIS 

In this section, we outline the quantities useful for de- 
tecting a possible floating KT phase and the adjacent KT 
phase transitions as well as the ordered phase. 
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A. Order parameters 

The complex XY order parameter deduced from the 
LGW analysis is given by: 



me* = (mi + m2e* ^ + mse* 



(8) 



where the m^ are the magnetisations of the three sub- 
lattices, and m is real and positive, m is close to 1 in 
the limit of zero temperature and vanishing transverse 
field. It equals zero for the disordered and KT phases, 
vanishing exponentially and algebraically in the limit the 
system size L — » cxd. There, the corresponding suscepti- 
bility is 

X = L^rr?)lT. (9) 
To detect clock symmetry breaking, we consider 
(m6cos(60)) 

It is easy to check that cg equals zero for disordered and 
KT phases, cg equals —1 for the (+0—) phase, and Cg 

equals 1 for the (H ) phase. We have chosen not 

directly to average cos(66') as its value fluctuates most 
strongly when m is small, that is to say, the ordering we 
are trying to determine the details of is weakest. The 
exponent m^ has been chosen, somewhat arbitrarily, as 
the corresponding term in the LGW action is -0+ + V'- = 
m^ cos(60). 



B. Binder cumulant 

The appropriate Binder cumulant isi^l 

U^l-J^. (11) 



3(m2)2 

The Binder cumulant has a scaling dimension of zero. 
It thus has the advantage of not requiring fitting an un- 
known leading exponent. In the limit L — > cx), the Binder 
cumulant has the following behaviour: Ul — > at disor- 
dered phase, Ul 2/3 at ordered phase, and Ul U* 
at a critical point. For an extended critical phase, the 
value of U* can drift, so that its value depends on the 
precise location within the critical phase. 

Since the Binder cumulant has a scaling dimension of 
zero, the curves for different system sizes at a critical cou- 
pling should fall on a line of points where Ul — Ul' (at 
least in the region where corrections to scaling are small) . 
As we describe below, this criterion in fact provides an 
overestimate of the size of the critical phase. 



C. Locating the KT transitions 

Here, we briefly describe the finite-size scaling analysis 
appropriate for KT transitions. This analysis follows that 



used by Challa and Landaui^^ In order to determine the 
presence of a KT phase, we check if we can fit our sim- 
ulations to the scaling forms predicted by KT theoryiSS 
This has the advantage of restricting the number of fit- 
ting parameters, which in the most general case (allowing 
for critical exponents different from the KT ones) would 
be too large to be practical. 

In particular, the correlation length ^ and the suscep- 
tibility X behave as 



^ (X exp(at "^/^), 
m cx r"^^ 



(12) 



where a is a nonuniversal constant, t = (T — Tc)/Tc is a 
reduced temperature, and is a critical temperature. 

The finite-size scaling form of order parameter and sus- 
ceptibility is given by 



mL = L ''mo(^/L), 
XL = L^Xoi^/L), 



(13) 



where mo and Xo are unknown universal functions and 
b and c are constants. It follows from (|12|l that in the 
infinite system size limit, we have b = rj/2 and c = 2 — rj. 
Therefore, at a critical point, one has: 



niL cx 
XL (X L'-\ 



(14) 
(15) 



If we have an extended critical phase, these relations 
should hold over finite temperature range from the upper 
critical temperature T2 to the lower critical temperature 
Ti. The critical exponent ry should vary continuously 
from T2 to Ti. Plotting In(mi) or \ii{xl) versus ln(L), 
we can find the critical exponent 77 at any point of the 
extended critical phase. 
We can rewrite Eq. H13|) as 



tulL'' = mo(L"^ exp(ai"^/^)), 
XLL-'' = Xo{L-'exp{at-'/')). 



(16) 
(17) 



Eq. (O is valid for T < Ti and ^ is valid for T > T2. 
For an appropriately chosen set of parameters a, c, Ti 
the plot of xlL~'^ versus exp(ai~^/2) should col- 
lapse onto a universal curve for different system sizes 
L. The same should hold for a plot of ttilL'' versus 
L^^ exp(af~^/^). From such a fit, the upper and lower 
critical temperatures can be determined. 



V. RESULTS 

In Fig. n we present the phase diagram that we find 
from Monte Carlo simulations. There are three phases: a 
disordered phase at high temperatures (T J) or large 
magnetic field strengths F ^ J, an extended KT phase at 
intermediate temperatures, and an (-1-0—) ordered phase 
at low temperatures. 
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FIG. 1: The phase diagram of the transverse field Ising model 
on the triangular lattice. A fioating KT phase separates the 
ordered clock phase from the paramagnet (PM) . QCP denotes 
the location of the quantum critical point. 



A. The KT phase at T/J = 0.4 

Our most complete data set was taken for F/J = 0.4, 
for which we describe our analysis in detail in the next 
few paragraphs. We have chosen this value of F as a com- 
promise between the following requirements. Firstly, we 
want to stay well clear of the zero temperature critical 
points at F = and larger F = Fc and concomitant pos- 
sible crossover phenomena. Secondly, a high transition 
temperature is needed to allow us to simulate systems 
with only moderate extent in the imaginary time direc- 
tion. Thirdly, as our algorithm slows down as the density 
of domain walls in the imaginary time direction increases, 
we would like to choose weak quantum fluctuations, that 
is to say, small values of F. V / J = 0.4, located left of cen- 
tre of the KT dome (Fig.^ thus appears to be a sensible 
choice. 

First, Fig. 121 shows the behaviour of the order param- 
eter m as a function of inverse temperature [3 for differ- 
ent system sizes. One can notice easily that the order 
parameter has strong dependence on the system size in 
wide range of temperatures - it decreases with increasing 
system size. 

Fig. O is a doubly logarithmic plot of the behaviour 
of the order parameter as a function of the system size. 
At high temperatures, it extrapolates to zero, whereas 
it levels out to a constant at low temperatures. For a 
wide range of temperatures in between, the data points 
display linear behaviour (without any evidence of log- 
arithmic corrections), beginning at (3 J = 2.8 down to 
(3 J = 9.0. The exponent rj varies continuously from 0.323 
at 15 J = 2.8 to 0.092 at /3J = 9.0. 

This is indeed close to the range expected from KT 
finite-size scaling (|15f) . although there is an overestimate 
of the size of the KT phase when compared to the ex- 
pected range of critical exponents, which is between 1/4 
and 1/9. As in the case of the Binder cumulant below, 
this is a simple consequence of the fact that the differ- 
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FIG. 2: The order parameter m versus /3J for different system 
sizes. 



ence between the ordered and the critical phase is indis- 
cernible when the size of the system is much less than the 
correlation length. As the latter grows near the transi- 
tion, increasingly large finite-size systems in the ordered 
phase appear to be critical. 




FIG. 3: The order parameter m versus system size L at vari- 
ous temperatures. The lines connect the points to guide the 
eye. 

To check whether this critical region is indeed delin- 
eated by KT transitions, in Fig. 0] we show the data 
collapse (Eq. I17|l for the upper transition. The data 
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scales very well with the following set of parameters 
(32J — 3.5, a — 5.36, c — 1.736. Thus we can conclude 
that the system has a transition between the disordered 
phase and KT phase at J — 3.5 ±0.2. The critical 
exponent 77 — 0.263 ± 0.015 at the transition point. This 
value of the critical exponent 77 is close to the theoretical 
prediction 1/4, and the overestimate appears to be part 
of a systematic trend discussed below. 
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FIG. 4: Data collapse of the susceptibility for the upper tran- 
sition. 



ture approaches zero. This implies the existence of the 
(+0— )-phase at low temperatures for F/J = 0.4. 

We can also determine the lower transition tempera- 
ture T\ from the criterion of 77 = 1/9 at the lower tran- 
sition point, which yields fi\J — 7.5 ± 0.5. At the tran- 
sition temperature thus determined, one finds a crossing 
of cos(66') as a function of /? for different system sizesi^ 
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FIG. 6: The order parameter cos(6^) versus /3J for different 
system sizes. 



In Fig. |31 we show the lower temperature scaled data 
(see Ea. ll6|) . The data scales quite well (but not perfectly, 
and over a narrower range than for the upper transition) 
with the following set of parameters: (i\J = 8.0, a = 
1.2, c = 0.105. The error in determining the critical tem- 
perature is larger than in the high temperature case. We 
can conclude that the system has a transition between 
the KT phase and the ordered phase at /?i J = 8.0 ± 1.0. 
The critical exponent r\ = 0.105 ± 0.02 at the transition 
point. This value of the critical exponent r\ is again close 
to the theoretical prediction 1/9. 

r/J=0.4 2b=0.105 a=1.2 |3;,J=8 
1.14 I , , 1 




1.04 I ' ' 1 

0.1 1 10 

L"^exp(at"^'^) 

FIG. 5: Data collapse of the order parameter m for the lower 
transition. 

To determine the nature of the ordered phase, we con- 
sider the behaviour of cos(66') (Eq. I10|l . The plot of 
cos(60) as a function of the inverse temperature /? is 
shown in Fig. El cos(60) goes to —1 as the tempera- 



Next, we consider the flow diagram of the Binder cu- 
mulant, depicted in Fig.|7| In the high (low) temperature 
phase, its value approaches (2/3) as L — s- 00, whereas 
it levels off to a value V* in the KT phase which depends 
on the strength of the coupling. Challa and Landau2& 
proposed to use this levelling off as a diagnostic for the 
KT phase. We do indeed observe such a levelling over a 
wide range of temperature at T / J = 0.4 J. 

This diagnostic again systematically overestimates the 
extent of the critical phase. This is evidenced by the 
uppermost curves in Fig. [7| There is an inflection point 
at large system size where the value of Ul, having ap- 
parently levelled off, starts increasing again. As one ap- 
proaches the transition from the ordered phase, this point 
of inflection wanders to increasingly larger system sizes 
and hence beyond the scope of the simulations. 

B. Other values of V 

We now describe what we have found for other values of 
F. At small F/J <C 1, T/ J <C 1, J drops out as an energy 
scale as the Hilbert space is restricted to the classical 
ground states, with excited states frozen out by an energy 
gap of 0{J). The only dimensionless parameter is then 
/3F, so that the lines emanating from the origin of the 
phase diagram are straight. We estimate their slopes to 
be T^rp/r = 0.41 ± 0.05 for the lower transition, and for 
the upper: TJ^^/r = 0.86 ± 0.09 

The higher temperature data collapse does not al- 
ways give the critical exponents 77 at critical tempera- 
tures close to 1/4, especially for large F. For example, 
we get the following upper set of critical temperatures 
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FIG. 7: The Binder cumulant Ul versus system size L at 
various temperatures. F/J = 0.4. Note the different scales 
for the ordinate. 



FIG. 8: The Binder cumulant Ul versus system size L at 
various temperatures. F/J = 0.8. Note the different scales 
for the ordinate. 



and critical exponents 77: (3cJ — 5.8(6), 77 — 0.287(25) at 
F/J = 0.2, f3cJ = 3.5(2), 77 = 0.263(15) at T/J = 0.4, 
PcJ = 2.6(2), 77 = 0.230(20) at F/J 0.8, and (3cJ = 
3.0(3), 77 = 0.203(25) at F/J = 1.2. For T/J = 0.2, the 
high-temperature data scales very well over a wide range 
of temperatures. Thus it is difficult to determine the 
precise value of the upper transition temperature at this 
F. 

The systematic trend of decreasing 77 with increasing 
F is in accordance with the fact that 77 at the zero tem- 
perature transition at large F is different. For the d — 3 
XY universality class, 77 is in fact close to so that the 
increasing proximity of this fixpoint should be expected 
to show up in a correction of this kind. One can account 
for the systematic trend of increasing 77 with decreasing F 
in the same way. Indeed, the frustrated triangular Ising 
model has a critical point at T = 0, where 77 — 1/2M 

Probably related to this crossover is the fact that, as F 
increases, the plots of the Binder cumulant as a function 
of system size fail to display the clear flattening visible 
in Fig. [7| for T/J — 0.4 until a system size which is sub- 
stantially larger; this is displayed for T / J — 0.8 in Fig.|Sl 

The lower transition at larger F poses a problem for the 
alogrithm we use. Since we used linked lists to store the 
location of the domain walls in the temporal direction, 
the algorithm slows down considerably as F increases. 
Together with the increasingly slowly decaying correla- 



tions in the space direction as the temperature is lowered, 
this leads to a considerable increase in the correlation 
times of our algorithm. As a result, the data we were 
able to obtain were too noisy to permit a useful data col- 
lapse. We can, nonetheless, try to identify the location 
of the transition by tracking where 77 crosses 1 /9 (Fig. ^ 

We note that the classical six-state clock model, which 
has a phase diagram very similar to the one we find here, 
has been studied by many authorsiS2iS As is the case for 
us, almost all the authors report that it is much more 
difhcult to get data collapse at the lower KT transition 
than at the higher. Only in a recent paper, Tomita 
and Okabe^^ using a 'probability-changing cluster al- 
gorithm', have found the lower transition easily. Given 
the frustrated and higher dimensional nature of our spin 
model, such a cluster algorithm is not available to us. 

The nature of the ordered phase at low F is of the 
(-1-0—) type. There is in principle the possibility of a 
change of sign of vq for entropic reasons as the couplings 
are variedAiS As F increases, it becomes increasingly 
hard to determine the sign of cos(66') for the system sizes 
available to us. The largest F for which we can con- 
fidently state that the lower KT transition is into the 
(-1-0—) phase is F/J = 1.2. Although this difficulty may 
be in part due to a decrease in the strength of V(i, there 
is no supporting numerical evidence that the sign of iig 
ever changes. 
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C. The quantum critical point 

We now turn to the properties of the quantum crit- 
ical point, i.e. to the zero temperature transition from 
the clock symmetry broken ordered phase to the para- 
magnet. We can determine its approximate location to 
be at Tf,/ J — 1.65 ± 0.05. For large F, the magnet is in 
a quantum paramagnetic phase with a size-independent 
susceptibility. 

We have not directly attempted to determine the prop- 
erties of this transition as extensive previous simula- 
tions on the classical stacked magnet exist. At the 
time, there were some suggestions that the observed crit- 
ical exponents were in fact not those of a d = 3 XY 
model^'* and perhaps altered due to the presence of an- 
other instabilityiiSii^ 

The structure of our phase diagram lends support to 
the d = 3 XY universality class scenario via an indepen- 
dent route. By inducing a dimensional crossover through 
switching on a finite temperature, we find that the highly 
nontrivial phase diagram is that predicted by the same 
Landau theory which gives the d = 3 XY universality 
class. This diagnostic is perhaps more robust than a di- 
rect determination of the critical exponents, which can 
be influenced by corrections to scaling or the proximity 
of other instabilities. 

The shape of the phase boundary near the quantum 
critical point follows from the knowledge of the criti- 
cal exponents of the quantum phase transition.— The 
boundaries of the KT phase near Fc follow the trajec- 
tory Tkt oc |F — Fcl"^^, where, for the present case, the 
dynamical critical exponent z = 1 and v is close to 2/3 

VI. CONCLUSION 

We have demonstrated that the common action of 
thermal and quantum fluctuations in the triangular lat- 
tice transverse field Ising model generate an interesting 
fluctuation-driven phase diagram including an extended 
critical phase bordered by a pair of KT transitions. We 
have employed several diagnostics for the presence of the 
KT phase and our results are consistent in considerable 
detail with what one would expect from an analysis based 
on a d = 2 + 1 dimensional XY symmetric action with 



^ For a non-technical introduction to frustrated magnets, 
see R. Moessner, Can. J. Phys. 79, 1283 (2001); re- 
views of exact diagonalisations and experiments, respec- 
tively, are C. Lhuillier, P. Sindzingre, and J.-B. Fouet, 
cond-mat/0009336 and P. Schiffer and A. P. Ramirez, 
Comments Cond. Mat. Phys. 18, 21 (1996). 

^ J. Villain, R. Bidaux, J. P. Carton, and R. J. Conte, J. 
Phys. Paris 41 1263 (1980). 

3 E. F. Shender, Sov. Phys. JETP 56, 178, (1982). 

* C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 



a six-fold clock anisotropy for a temporal dimension of 
tunable extent. This provides indirect support for the 
XY nature of the quantum critical point. The dynami- 
cal properties of the critical phase in close proximity to 
the quantum critical point are a subject worthy of a study 
in their own right. 

The agreement between numerics and theory, however, 
is not complete. There are some deviations from the pre- 
dicted critical exponents, which are probably due to a 
crossover to the zero temperature critical points. As we 
have used several diagnostics to determine the nature 
and location of phases and transitions, we are in a posi- 
tion to compare their relative reliability and find that, for 
the phase transitions in particular, considering the sys- 
tem size independence of the Binder cumulant and the 
power-law dependence of the magnetisation on system 
size systematically overestimate the extent of the critical 
phase. 

In summary, the results we have obtained demonstrate 
once again that frustrated magnets provide a good start- 
ing point for finding unconventional phases and phase di- 
agrams. In this particular case, by using a tunable com- 
bination of thermal and quantum fluctuations, we have 
managed to realise a standard model from statistical me- 
chanics, the XY model with sixfold clock anisotropyjiS 
in dimensions d = 2 and d = 3, in terms of another 
one, namely the Ising model on the triangular lattice. 
Given other magnets in this class realise unusual order 
parameters,* this approach should provide more opportu- 
nities for studying exotic phase diagrams based on simple 
model spin systems. 



Acknowledgements 

We would like to thank D. Huse, A. Karlhede and es- 
pecially S. Sondhi for many useful suggestions and dis- 
cussions throughout the course of this work. RM is also 
grateful to P. Chandra and S. Sondhi for collaboration on 
closely related work. In addition, we thank A. Bunker, 
S. Fujiki, T. H. Hansson, J. Lidmar, and M. Oshikawa 
for useful discussions. This work was in part supported 
by the Ministere de la Recherche et des Nouvelles Tech- 
nologies with an ACI grant. 



^ see, e.g. K. Damle and S. Sachdcv Phys. Rev. B 56, 8714 
(1997) and S. Sachdev, cond- mat /9508080 for a pedagogi- 
cal review. 

® J. Brooke, T.F. Rosenbaum, and G. Aepph, Nature 413, 

610 (2001); G. Aeppli, talk at conference on the Physics of 

Frustration, Santa Fe, 2002 (unpublished). 

G. H. Wannier, Phys. Rev. 79, 357 (1950); R. M. F. 

Houtappel, Physica 16, 425 (1950). 
* R. Moessner, S. L. Sondhi, and P. Chandra, Phys. Rev. 

Lett. 84, 4457 (2000), Phys. Rev. B 64, 144416 (2001); 



9 



R. Moessner and S. L. Sondhi, Phys. Rev. B 63, 224401 
(2001). 

^ D. Blankschtein, M. Ma, A.N. Berker, G.S. Grest and CM. 

Soukoulis, Phys. Rev. B 29, 5250 (1984). 
^° J. V. Jose, L. P. Kadanoff, S. Kirkpatrick and D. R. Nelson, 

Phys. Rev. B 16, 1217 (1977). 
" P. Chandra, P. Coleman and L. B. loffe, Phys. Rev. B 49, 

12897 (1994). 

^2 O. Heinonen and R. G. Petschek, Phys. Rev. B 40, 9052 
(1989). 

R. R. Netz and A. N. Berker Phys. Rev. Lett. 66, 377 
(1991). 

Alex Bunker, B. D. Gaulin, and C. Kallin, Phys. Rev. B 
48, 15861 (1993). 

M. L. Plumer and A. Mailhot, Physica A222, 437 (1995). 
For a review, see M. F. Collins and O. A. Petrenko, Can. 
J. Phys. 75, 605 (1997). 
" M. Suzuki, Prog. Theor. Phys. 46, 1337, (1971); ibid 56, 
2454 (1976). 

1* H. F. Trotter, Proc. Am. Math. Soc. 10, 545 (1959). 
^® N. Todoroki, Y. Ueno, and S. Miyashita, Phys. Rev. B 66, 
214405 (2002). 

D. J. Amit, Field Theory, the Remormalization group and 

Critical Phenomena, World Scientific (1984). 

B. B. Beard and U.-J. Wiese, Phys. Rev. Lett. 77, 5130 

(1996). 

^2 N. V. Prokof'ev, B. V. Svistunov, and I. S. Tupitsyn, 



Pis'ma Zh. Eks. Teor. Fiz. 64, 853 (1996) [JETP Lett. 64, 
911 (1996)]; Zh. Eks. Teor. Fiz. 114, 570 (1998) [JETP 
87, 311 (1998)]. 

H. Rieger and N. Kawashima, Eur. Phys. J. B 9, 233, 
(1999). 

^'^ R. H. Swendsen and J.-S. Wang, Phys. Rev. Lett. 58, 86 
(1987). 

U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 
^® For a pedagogical introduction, see W. Krauth, 
cond-mat/9612186 

K. Binder, Z. Phys. B 43, 119 (1981). 
^® M. S. S. Challa and D. P. L. Landau Phys. Rev. B 33, 

437 (1986). 
23 J. Kosterlitz, J. Phys. C 7, 1046 (1974). 
^° Indeed, M. Oshikawa (Phys. Rev. B 61, 3430 (2000)) notes 

that cos(6S) has a scaling dimension of zero for the six-state 

clock model. 

R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 
^2 J. Stephenson, J. Math. Phys. 11, 413 (1970). 

S. Fujiki, K. Shuton, S. Inawashiro, Y. Abe, and S. Kat- 
sura, J. Phys. Soc. Jpn. 55, 3326 (1986); A. Yamagata 
and I. Ono, J. Phys. A 24, 265 (1991); S. Fujiki and T. 
Horiguchi, J. Phys. Soc. Jpn. 64, 1293 (1995). 
Y. Tomita and Y. Okabe, Phys. Rev. B 65, 184405 (2002). 
For reviews, see S. L. Sondhi, S. M. Girvin, J. P. Carini 
and D. Shahar, Rev. Mod. Phys. 69, 315 (1997), or Ref.H 



